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ABSTRACT: This is a numerical analysis of, forced convection heat transfer inside a channel of rectangular 
section, containing some inclined baffle plates. We have used a numerical model based on a finite-volume 
method, and The SIMPLE algorithm [6] was adopted to solve the coupling pressure-velocity problem. Effects of 
various parameters of the baffles, such as, baffle 's inclination angle for two cases 'channel with one or two 
baffles ' on isotherms, streamlines, temperatures distributions and Nusselt number values was studied. It is 
concluded that baffles angle inclination has a meaningful effect on isotherms, streamlines and total heat 
transfer through the channel. Indeed, the increase in the baffle inclination improves heat transfer within the 
channel. Also, heat transfer becomes decreasingly important with adding baffles. 
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I. INTRODUCTION: 

Turbulent forced convection in heat exchangers with different type of baffles was investigated, 
recently, in a large number of experimental and numerical works [1-4]. This interest is due to the various 
industrial applications of this type of configuration such as cooling of nuclear power plants and aircraft engine 
... etc. S. . V Patankar and EM Sparrow [1] solved numerically the problem of fluid flow and heat transfer in 
fully developed heat exchangers. These one was equipped by isothermal plate placed transversely to the 
direction of flow. They found solid plates caused strong recirculation zones in the flow field. They concluded 
that the Nusselt number depends strongly on the Reynolds number, and it is higher in the case of fully 
developed then that of laminar flow regime. 

Demartini et al [2] conducted numerical and experimental investigations on turbulent flow inside a 
rectangular channel containing two rectangular baffles. They found that numerical results were in good 
agreement with those obtained by experiment. In conclusion, baffles play an important role in the dynamic 
exchangers studied. Indeed, regions of high pressure 'recirculation regions' are formed nearly to chicanes. 
Recently, Nasiruddin and Siddiqui [3] studied numerically effects of baffles on forced convection flow in a heat 
exchanger. The effects of size and inclination angle of baffles were detailed. They considered three different 
arrangements of baffles. They found that increasing the size of the vertical baffle substantially improves the 
Nusselt number. However, the pressure loss is also important. For the case of inclined baffles, they found that 
the Nusselt number is maximum for angles of inclination directed downstream of the baffle, with a minimum of 
pressure loss. 

More recently, Saim et al [4] presented a numerical study of the dynamic behavior of turbulent air flow 
in horizontal channel with transverse baffles. They adapted numerical finite volume method based on the 
SIMPLE algorithm and chose , k-s model for treatment of turbulence. Results obtained for a case of such type, 
at low Reynolds number, were presented in terms of velocity and temperature fields. They found the existence 
of relatively strong recirculation zones near the baffles. The eddy zones are responsible of local variations in the 
Nusselts numbers along the baffles and walls. 

We know that the primary heat exchanger goal is to efficiently transfer heat from one fluid to another 
separated, in most practical cases, by solid wall. To increase heat transfer, several approaches have been 
proposed. We can cite the specific treatment of solid separation surface (roughness, tube winding, vibration, 
etc.). This transfer can also be improved by creation of longitudinal vortices in the channel. These eddies are 
produced by introducing one or more transverse barriers (baffle plates) inside the channel. The formation of 
these vortices downstream of baffles causes recirculation zones capable of rapid and efficient heat transfer 
between solid walls and fluid flow. It is this approach that we will follow in this study. Indeed, we are interested 
in this work on the numerical modeling of dynamic and thermal behavior of turbulent forced convection in 

12 



Numerical Analysis of heat transfer in a channel wit incl. 



horizontal channel where two walls are raised to a high temperature. This channel may contain one or several 
rectangular inclined baffles. 

A special interest is given to the influence of different parameters, such as number and baffle inclination on 
heat transfer and fluid flow. 



II. MATHEMATICAL FORMULATION 

The geometry of the problem is shown schematically in Figure Lit is a rectangular duct with 
isothermal horizontal walls, crossed by a stationary turbulent flow. The physical properties are considered to be 
constants. 




At each point of the flow the velocity has components (u, v) in the x and y directions, and the 
temperature is denoted T. The turbulence modeling is handled by the classical model (k-e). k is the turbulent 
kinetic energy and 8 the viscous dissipation of turbulence. The Transport equations (continuity, momentum, 
temperature, turbulent kinetic energy and dissipation of turbulence) governing the system, are written in the 
following general form : 
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Where p is the density of the fluid passing through the channel and <|), Fcp and Sep are given by: 
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With : 
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u and ut represent, respectively, the dynamic and turbulent viscosities. 
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The constants used in the turbulence model (k-£) Are those adopted by Chieng and Launder (1980) [8] 
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Boundary conditions : 
At the channel inlet: 

u = U . ; v = 0 and T = T in 

k = 0,005 U * n and e = 0,1k 3 ' 2 
At solid walls: 

u = v = o .k = s = 0 and T = T w >■ T i„ 
At the channel exits: 

The gradient of any quantity, with respect to the longitudinal direction x is nul. Our goal is to determine velocity 
and temperature fields, as well as the turbulence parameters. Particular attention is given to the quantification 
parameters reflecting heat exchanges such as the Nusselt number (local and average). 

III. NUMERICAL PROCEDURE 

The computer code that we have developed is based on a finite volume method. Computational domain 
is divided into a number of stitches. To choose the number of cells used in this study, we performed several 
simulations on a channel (with or without baffles). Finally, we have opted for a 210 x 90 meshes. 
The mesh dimensions are variables; they tighten at the solid walls neighborhoods. Consequently, the stitch 
density is higher near the hot walls and baffles. 

We consider a mesh having dimensions Ax and Ay. In the middle of each volume we consider the 
points P, called centers of control volumes. E, W, N, S are the centers of the adjacent control volumes. We also 
consider centers, EE; WW, NN, SS. The faces of each control volume are denoted e, w, n, s 'figure 2'. 
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Figure 2: Mesh with P at center 
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By integrating the transport equation (1) of § on the control volume, we find, a relation to the x direction: 

\ d <j>^\ \ d </> 1 rl 
L/j u </> \ e - [p u 0 \ w = yT — — j - yT —— j +[_S </I \ p 

Where [s ^ ] is the source term. 

One note F = pU, the convection flow and D = T I A x , the diffusion coefficient. And with taking : 

F = \p u 1 ; F = \p u 1 

D =[r/Axl ;D =[r/Axl 

And: = </> E ; [d^] w = </> P - K 

We found : F ^ e - F w <P w = D e {<f> E - </> p ) - D w {</> p - <f, w ) + [s J p 

To estimate $ on the faces " e " and " w " we opted for the classical quick scheme [5] which is a quadratic 
forms using three nodes. The choice of these nodes is dictated by the direction of flow on these 
faces (u >- 0 or u -< 0) . 

Finally, the transport equation (1) is discretized on the mesh with P at center as : 

a p </> p =a w </> w +a E <f> E + a ww </> ww + a EE 4> EE + [s , ] p (2) 

f 6 13 
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With : 
Where 



a = 1 or 0 , if F >■ 0 or F -< 0 
a = 1 or 0 , if F ^0«r F ^0 



The same thing is done for the vertical "y" direction, using the faces "n" and "s" and by introducing the 
Quick diagram nodes N, NN, S and SS. 

If convection flow and diffusion coefficient of the sides "e", "w", "s", "n" are known And especially 

the source term [s ] , the solution of equation (2) then gives us the values of § in different P nodes. 

One notes that to accelerate the convergence of equation (2) we have introduced a relaxation factor : 

a P r i 1 - P o 

— <t> P = a w (f> w + a E (/> E + a ww 0 WW + a EE </> EE + [S^\ p + - a p (/> p (3) 

Where ^° is the value of </> p in the previous step. 

The source term appears especially in the conservation of momentum equations in the form of a 
pressure gradient which is in principle unknown. To get around this coupling we have chosen to use in our code 
the "SIMPLE" algorithm developed by Pantankar [7]. The basic idea of this algorithm is to assume a field of 
initial pressure and inject it into the equations of conservation of momentum. Then we solve the system to find a 
field of intermediate speed (which is not fair because the pressure isn't.) The continuity equation is transformed 
into a pressure correction equation. This last is determined to find a pressure correction that will inject a new 
pressure in the equations of motion. The cycle is repeated as many times as necessary until a pressure correction 
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equal "zero" corresponding to the algorithm convergence. In the end we solve the transport equations of T, k 
and £. In this approach, a problem is encountered. It is known as the checkerboard problem. 

The risk is that a pressure field can be highly disturbed by the sensed formulation which comprises 
performing a linear interpolation for estimate the pressure value on the facets of the control volume. To 
circumvent this problem we use the so-called staggered grids proposed by Harlow and Welch [6]. In this 
technique, a first grid pressure (and other scalar quantities T, k and £) is placed in the center of the control 
volume. While other staggered grids are adopted for the velocity components u and v (see Figure 3). 
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The scalar variables, including pressure, are stored at the nodes (I, J). Each node (I, J) is surrounded by 
nodes E, W, S and N. The horizontal velocity component u is stored on faces "e" and "w", while the vertical 
component v is stored on the faces denoted "n" and "s". 

So the control volume for pressure and other scalar quantities T, k and £ is 

[(« + 1, j); (i + l, ; + l); (i, j + 1); (i, j)] 

For component u control volume is 
[(/, ;); (/, j + 1); (/ - 1, j + 1); (/ - 1, j)] 

While for the v component we use: [(;' + 1, J - 1); (i + 1, /); (i, J); (i - 1, / - 1)]. 

By integrating the conservation momentum equation in the horizontal direction on the 
volume[(/, ;'); (/, j + 1); (/ - 1); (/ - 1, ./)], we found : 

a i,j u i,j = 2 a n b u n b + l> ( 7 - !. J ) - p ( 7 - J )]( y y +J - y ) ( 4 ) 

Where: ^ a nb u nb = a ww u ww + a w u w + a EE u EE + a E u E + a ss u ss + a s u s + a m u m + a N u N 

The coefficients anb are determined by the quick scheme mentioned above. 

Similarly, the integration of conservation momentum equation in the vertical direction on the volume 

[(» + 1,7- 1); (i + 1,J); (i,J ); (i - 1,J - l)], gives: 
a .,i v .,i = X a n», v n b + [PV.J ~ 1) - P(I,J)](x M - x,) (5) 

One considers primarily an initial pressure field P*. The provisional solution of the equations (4) and (5) will be 
denoted u* and v*. We note that u* and v* doesn't checks the continuity equation, and: 

a^u'xj = X a nl u' nb + [p\l - 1,7) - /»*(/, 7 )](y j+l - y .) (6-a) 

a = X a „ b v "" b + [p'V-J -D- - *,) (6-b) 

At this stage any one of the three variables is correct, they require correction : 
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\ p'= p - p" 
I 

■j u ' = u - u (7 ) 

v'= v - V 

Injecting (7) in equations (4), (5) and (6) we find : 

«,.,»'../ = £ a n„ u 'n b + i p '( I -l,J)~ P'(/,/)](y y+1 - y,) (8-a) 
«i.y v '/.; = Z a rt v' rt +[P'(/, / -l)-P'(/,/)]U (+1 - *,) (8-b) 

At this level, an approximation is introduced. In order to linearize the equations (8), the terms 

Za .»' , V <j . v' 

are simply neglected. Normally these terms must cancel at the procedure convergence. That is to say that this 
omission does not affect the final result. However, the convergence rate is changed by this simplification. It 
turns out that the correction P' is overestimated by the Simple algorithm and the calculation tends to 
diverge. The remedy to stabilize the calculations is to use a relaxation factor. 

We Note that a further treatment of these terms is proposed in the so-called algorithms "SIMPLER" and 
"SIMPLEC" [7]. 



The equations (8) become : 

u \ } = d u [P'(I - 1,7) - P'(I,J)] (9 - a) 

v', . = d 1} [P'(I,J - 1) - P'(/,7)] (9 - b) 



With : a " and 



Equations (9) give the corrections to apply on velocities through the formulas (7). We have therefore : 

u , = u'i.j + d., \P'(I - 1,7) - P '(7,7)1 

J (10-a) 

v, . = v'i.j + d, \P'(I,J - 1) - P'(/,7)l 

J (10-b) 

Now the discretized continuity equation on the control volume of scalar quantities is written as 
follows : 

[{puA) MJ - (p«A).J- [(pvA) Ij+l - (pvA), .]= 0 
Where A are the size of the corresponding faces. 

The introduction of the velocity correction equations (10) in the continuity equation gives a final equation 
allowing us to determine the scope of pressure corrections P': 

a' 7/ P'(/,7) = a' [+lJ P'(I + 1,7)+ a',_ u P'{I -1,7)+ a' IJ+1 P'(I,J + l) + a' IJ _ J P'(/,7 -1) (11) 

The algorithm can be summarized as follows : one starts from an initial field p * , u * , v * and </> * , with ^ 

represents the scalar quantities T, k and £. Then the system (6) is solved to have new values of u "and v * . 
Then the system (1 1) is solved for the corrections field pressure P. 

Thereafter, pressure and velocity are corrected by equations (7) and (10) to have P, u and v. 
We solve the following transport equation of scalar quantities (2), § — T, K and £ . 

Finally, we consider: P* = P ; u* = u ; v * = v ; ^ > ' = (j> . 

and the cycle is repeated until convergence. 

IV. RESULTS AND DISCUSSIONS 

The dimensions of the channel presented in this work are based on experimental data published by 
Demartini et al [2]. The air flow is carried out under the following conditions. 

> channel length: L = 0.554 m; > channel diameter D = 0.146 m; 

> baffle height : h = 0.1 m; > Baffle thickness: 5 = 0.01 m ; 

> Reynolds number: Re = 8.73 1 0 4 ; 
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The hydrodynamic and thermal boundary conditions are given by : 
At the channel inlet: 

u = u tn = 7.8m I s ; V = 0; T = T in = 300 °K 
On the channel walls: 
u = v = 0 and T = T w = 373 °K 

At the channel exit the system is assumed to be fully developed: 

8u / _ Qv / _ 8 T / _ q 
/ 8x / dx / Sx 

First, we compared the structure of streamlines and isotherms in channels with baffles having different 
inclination angles. Baffle is located at the abscise Lj = 0.218 m. Indeed, in figures 4 and 5 we present 
streamlines and isotherms for three values of inclination angle. 

These results clearly show the importance of presence of baffle (acting as a cooling fin) [9]. Since the 
recirculation zone, downstream of the baffle, depend on the inclination angle so this last affect isotherms and 
streamlines structures and affect, consequently, the heat transfer through the channel. Indeed the baffle increases 
heat transfer between the wall and the fluid. So this increase depends on baffle position [9] and angle 
inclination. 




Case 1-1 a = 45° 




(b) Case 1-2 a = 60° 




(c)Case 1-3 a = 90° 



Figure 4: streamlines 









Case 1-1 a = 45° 









Case 1-2 a = 60' 
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Case 1-3 a = 90° 



figure 5: Isotherms 'one baffle' 



It is clear that the isotherms are more condensed near baffle, top and bottom walls when the inclination 
angle a increases. This indicates an increase of heat transfer near baffle. Indeed, the increase in a expanded 
exchange area between fluid and walls. This is because the flow is blocked by the baffle when a is small. In 
addition, when a increases, the recirculation zone becomes increasingly important (see Figure 4). This causes an 
acceleration of turbulent flow in this region, which improves heat transfer within the channel. 
We present figures 6 and 7 in order to identify the influence of a on velocities and temperature profiles along y 
at x = 0.45 m ; from channel inlet. 
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Figure 6: Horizontal velocity Profiles along y at x = 0.45 m « influence of a » 

These results allow us to conclude that the effect of baffle inclination on temperature profile and 
velocity distribution depend on channel regions 'bottom and top'. Indeed, the velocity is maximal in the top half 
channel when a is maximal and the inverse is true in the bottom half channel. 
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Figure 7: Temperature profiles along y at x = 0.45 m « influence of a » 
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Concerning temperature profiles at x = 0.45 m, we can see that there are similar between cases 1-1 and 
1-2, with a nearly constant difference. However, the temperature is maximal for cases 1-1 and minimal for cases 
1-2. 

In the purpose of measuring the influence of ' a ' on local heat transfer, we have presented in figure 8, 
the local Nusselt number along the channel for the considered cases. 

350 i 1 




It shows that upstream of the baffle (x <0.15 m) curves are confused , while, just after this location, 
effect of baffle inclination on Nusselt number become increasingly important away from baffles. Indeed, the 
local Nusselt number in the case 1-1 is more important than the other cases. This indicates that heat transfer is 
higher in the case 1 - 1 . 




Case 2-1: a = 45° 




Case 2-2 : a = 60 




Case 2-3 : a = 90° 



Figure 9 : Streamlines 'two baffles' 

We consider two baffles of same height h = 0.1 m. The first baffle is placed at the distance L \ = 0.15 m, 
while the second is positioned at the distance d i = 0.10 m from the first. In order to investigate the inclination 
angle we have chosen three different inclination angles; a = 45° (case 2-1), a = 60° (case 2 -2) and a = 90° 
(case 2 -3). We note the existence of two recirculation zones downstream of the first baffle. Also the first 
recirculation zone delimited by baffles becomes increasingly important with increasing a, which contribute to an 
increase in heat transfer in this area, as shown in Figure 10. Indeed, when the baffles are vertical a = 90°, the 
fluid is blocked in the chimney delimited by two baffles. Which reduces the flow speed at that location and thus 
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there will be a decrease in the heat transfer in this area. In addition, when a increase, the fluid has sufficient 
space to move rapidly, hence heat transfer increases in this zone 'see figurel3 '. 












\ 






-v. ( 
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Case 2-3 : a = 90° 



Figure 10 : Isotherms 'two baffles' 

In the following figures 11 and 12, we compare the profile of horizontal velocity and temperature 
distribution. Calculation is done on the y-axis at x = 0.45 m from the channel inlet for three different baffles 
inclination angles. We can distingue two zones, the first defined by 0.02<Y<0.08 and the second by 
0.08<Y<0.02. In the first zone, we note that the increase of a contribute to an increase in the horizontal flow 
velocity; whereas, there is an opposite effect for the second region. 



Concerning temperature distribution at x = 0.45 m, one distingue two different zones: in the top region, 
corresponding to the interval 0.02<Y<0.08, fluid temperature is higher in the case of vertical baffles 'a = 90°'. 
Whereas it's lower in the bottom region, corresponding to -0.08<Y<0. This can be explained by the fact that the 
fluid passing the region between the first baffle and the channel top wall is heated more and more when the 
inclination angle, a, increases. This is because this region becomes small as the inclination angle a increases. 
Thereby the fluid reaches the second baffle with a higher temperature and velocity. Indeed the higher velocity 
allows the fluid to keep its temperature, due to the exchange with the lower cold region minimal. Consequently 
the fluid in the lower region has a minimal temperature. 
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Figure 11: Profiles of the horizontal velocity at x = 0.45 m « influence of a : 
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0.06 



3 0 2 3 0 4 3 0 6 3 0 8 3 1 0 3 12 3 14 3 1 G f (°K) 
Figure 12: Profiles of the total temperature at x = 0.45 m « influence of a » 
Figure 13 present the local Nusselt number along the channel for three considered inclination angles. 
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Figure 13: Local Nusselt number along the channel « Influence of a » 

Since the local Nusselt number shows the local heat transfer by convection, because it has a relation 
with the velocity and temperature distribution, its structure depend on the latter. Indeed, we distingue different 
zones. The first defined by 0 <x <0.35 m while the second by x>0.35 m. For the first region, the local nusselt 
number is higher for a = 60° and lower for a = 45° which means that the heat transfer is important in the first 
case and lower in the second case; whereas the reverse is true for the second zone. One note that the heat 
transfer is maximal for a = 45° and minimal for a = 60°. Let's note that local nusselt number for every case 
present two 'peaks' corresponding to baffles positions. 

We chooses to show the influence of baffles number on heat transfer along the channel for vertical 
baffles 'inclination angle a = 90°'. We have chosen to show case 1-3 'channel having single baffle' and case 2-3 
'channel containing two baffles', see figure 14. It should be noted that generally, heat transfer is inversely 
proportional to baffles number. Indeed, the Nusselt number characterizing heat transfer within the channel 
decrease with increasing baffles number. 
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V. CONCLUSION 

The thermal behavior of a stationary turbulent forced convection flow within a baffled channel was 
analyzed. The results show the ability of our code to predict dynamic and thermal fields in various geometric 
situations. We studied mainly the influence of baffles inclination and number on heat transfer and fluid flow. 
One can conclude that the increase in the baffle inclination improves heat transfer between channel walls and 
fluid passing through it; and that heat transfer becomes decreasingly important with adding baffles. 

In perspective, we intend deepen and clarify our results. Indeed, we will adapt our code to others 
geometric cases (non -rectangular baffles or baffles). Finally, we will also try to refine more the turbulence 
model. 
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